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Abstract 

In  this  study,  a  previously  developed  nickel-metal  hydride  (Ni-MH)  battery  model  is  applied  in  conjunction  with  experimental 
characterization.  Important  geometric  parameters,  including  the  active  surface  area  and  micro-diffusion  length  for  both  electrodes,  are 
measured  and  incorporated  in  the  model.  The  kinetic  parameters  of  the  oxygen  evolution  reaction  are  also  characterized  using  constant 
potential  experiments.  Two  separate  equilibrium  equations  for  the  Ni  electrode,  one  for  charge  and  the  other  for  discharge,  are  determined  to 
provide  a  better  description  of  the  electrode  hysteresis  effect,  and  their  use  results  in  better  agreement  of  simulation  results  with  experimental 
data  on  both  charge  and  discharge.  The  Ni  electrode  kinetic  parameters  are  re-calibrated  for  the  battery  studied.  The  Ni-MH  cell  model 
coupled  with  the  updated  electrochemical  properties  is  then  used  to  simulate  a  wide  range  of  experimental  discharge  and  charge  curves  with 
satisfactory  agreement.  The  experimentally  validated  model  is  used  to  predict  and  compare  various  charge  algorithms  so  as  to  provide 
guidelines  for  application-specific  optimization. 
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1.  Introduction 

The  nickel-metal  hydride  (Ni-MH)  battery  is  an  attrac¬ 
tive  replacement  for  the  Ni-Cd  battery  due  to  its  high  energy 
density,  high  rate  capability,  long  cycle  life  and  lack  of 
poisonous  materials.  The  demand  for  nickel-metal  hydride 
batteries  has  recently  grown  for  applications  ranging  from 
portable  electrics  to  electric  and  hybrid-electric  vehicles.  For 
clean  transportation,  Ni-MH  battery  is  presently  the  most 
promising  battery  for  hybrid-electric  vehicles  based  on  the 
performance  and  cost.  Extensive  efforts  have  been  made  to 
develop  various  advanced  nickel-metal  hydride  batteries  to 
meet  commercial  requirements. 

Traditionally,  experimental  testing  is  a  main  tool  to  test 
and  design  batteries.  However,  experiments  are  time  con¬ 
suming  and  costly,  and  experimentally  it  is  difficult  to 
determine  the  internal  process  during  charge  and  discharge. 
Modeling  and  simulation  of  batteries  is  a  useful  method  for 
battery  designers  because  a  good  cell  model  can  be  used  to 
identify  battery  mechanisms,  predict  the  cell  performance 
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for  design  and  optimization,  and  reduce  the  cost  and  time 
expended  on  the  experiments. 

Several  models  have  been  developed  to  study  the  nickel 
electrode  [1-5,13],  metal  hydride  electrode  [6-9]  as  well  as 
the  complete  nickel-metal  hydride  battery  [10,12,13].  Vidts 
and  White  [4]  developed  a  pseudo-2D  mathematical  model 
of  a  sealed  nickel-cadmium  cell  including  proton  diffusion 
and  ohmic  drop  in  the  nickel  electrode.  In  addition,  the 
oxygen  evolution  reaction  was  coupled  with  main  electro¬ 
chemical  reactions  on  both  electrodes.  A  sensitivity  analysis 
was  performed  to  determine  the  effects  of  various  para¬ 
meters  in  the  model.  They  found  that  at  the  beginning  of 
discharge,  the  sensitivity  of  exchange  current  density  of  the 
reduction  of  NiO(OH)  is  the  largest;  as  the  discharge 
proceeds,  the  sensitivity  coefficient  of  the  transfer  coeffi¬ 
cient  of  the  reduction  of  NiO(OH)  becomes  the  largest 
gradually;  the  effect  of  diffusion  coefficient  of  proton  on 
the  discharge  efficiency  is  the  largest  at  higher  discharge 
rates.  In  their  work,  however  no  experimental  data  was  used 
to  verify  the  model. 

Paxton  and  Newman  [10]  developed  a  nickel-metal 
hydride  cell  model  using  new  data  on  the  diffusion  coeffi¬ 
cient  of  proton  in  the  nickel  active  materials  as  measured 
by  Motupally  et  al.  [11].  In  this  model,  the  open-circuit 
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potentials  for  both  electrodes  were  experimentally  measured 
as  functions  of  the  state  of  charge.  They  argued  that  the 
proton  diffusion  could  be  neglected  in  the  low  rate  range  they 
studied.  They  showed  that  during  discharge,  the  potential 
drop  at  the  nickel  electrode  is  the  largest  loss  of  the  nickel- 
metal  hydride  cell.  Some  experimental  discharge  curves 
from  three  different  cells  running  at  the  0.5  C  discharge 
were  compared  with  their  model,  but  poor  agreement  can  be 
seen  due  to  a  lack  of  most  physico-chemical  properties  for 
the  cells.  The  oxygen  evolution  reaction  was  not  included. 

A  three-phase  electrochemical  model  was  developed  by 
Gu  et  al.  [12]  for  nickel-metal  hydride  cells  using  the  micro- 
macroscopic  coupled  approach.  By  including  the  oxygen 
reaction  in  the  cell,  the  effect  of  oxygen  evolution  on  the  cell 
voltage  was  assessed,  in  particular,  the  charge  acceptance  is 
correlated  with  the  charge  rate  and  the  evolution  of  oxygen  is 
responsible  for  the  cell  pressure  build-up.  In  addition,  they 
compared  experimental  data  and  modeling  results  for  dif¬ 
ferent  rate  charge  and  discharge  processes;  however,  the 
agreement  was  not  satisfactory  in  general  for  various  charge 
and  discharge  rates. 

The  objective  of  this  research  is  three-fold:  (1)  to  produce 
benchmark  experimental  data  with  all  needed  physico-che¬ 
mical  properties  that  will  be  useful  for  all  model  validation 
exercises;  (2)  to  refine  the  first-principles  mathematical  model 
of  Gu  et  al.  [12]  to  more  accurately  describe  the  behavior  of 
commercial  Ni-MH  cells  under  a  wide  range  of  operating 
conditions,  particularly  constant  current  (CC)  charge  and 
discharge,  and  constant  voltage  (CV)  operation;  and  (3)  to 
provide  some  design  change  examples  for  optimizing  the 
performance  of  the  battery  using  the  validated  model. 

2.  Model  development 

The  nickel-metal  hydride  battery  consists  a  Ni  positive 
electrode  and  a  MH  negative  electrode.  During  the  discharge 
and  charge  processes,  the  electrochemical  reactions  taking 
place  at  both  electrodes  are  given  by: 

Negative  electrode 


discharge 

MH  +  OH  «=>  H20  +  e~ 

charge 

(1) 

with  the  side  reaction 

5O2  +  H20  +  2e~  ->  20H- 

(2) 

H2  +  20H  2H20  +  2e 

(3) 

Positive  electrode 

discharge 

NiOOH  +  H20  +  e~  Ni(OH).,  +  OH 

charge 

(4) 

with  the  side  reaction 

OH  — >i02  +  H20  +  e- 

(5) 

2H20  +  2e~  ->  H2  +  20H 

(6) 

For  these  electrochemical  reactions,  we  follow  the  deriva¬ 
tion  of  Vidts  and  White  [4]  and  Vidts  et  al.  [6].  Both  liquid 
and  solid  phases  are  involved  in  species  transfer  during 
charge  and  discharge  processes,  and  O2  transport  in  the  gas 
and  liquid  phases  is  considered.  Charge  conservation  equa¬ 
tions  are  applied  to  both  solid  and  liquid  phases.  From  these 
conservation  equations  the  solid  and  liquid  phase  potentials 
can  be  determined.  The  model  used  in  this  paper  is  identical 
to  that  of  Gu  et  al.  [12]  and  hence  is  not  repeated  herein. 

A  battery  module  as  schematically  shown  in  Fig.  1  is 
composed  of  many  cells  connected  in  series,  as  well  as  in 
parallel.  The  module  modeling  can  be  approached  by  con¬ 
sidering  cells  to  be  connected  in  series  and  in  parallel  as  well 
as  resistances  of  inter-cell  and  terminals.  The  battery  module 
voltage  can  thus  be  related  to  the  cell  voltage  by  the  following: 

^  ~  E^conn  “  E^t  (7) 

7=  1  7=  1  7=1 

where  Ns  is  the  number  of  cell  in  series,  Is  is  the  current 
flowing  through  the  battery  stack,  RJconn  is  the  j-th  connector 
resistance  and  R[  is  the  j-th  terminal  resistance.  Assuming 
that  all  cells  in  the  module  behave  identically,  the  cell 
voltage,  V{  =  14,  can  thus  be  expressed  as: 

Vc  =  ^  [Vs  +  (Ns  -  1  )RconnIs  +  2 RJS]  (8) 

P's 

Here,  14  is  the  common  voltage  of  single  cells. 

3.  Numerical  procedure 

As  elaborated  in  Gu  et  al.  [12],  all  equations  in  this  model 
are  discretized  by  the  finite  volume  method  and  solved  by  a 
general-purpose  computation  fluid  dynamic  (CFD)  code. 
The  initial  condition  is  set  according  to  the  initial  state  of 
charge  (SOC)  and  species  concentrations.  The  boundary 
condition  depends  on  the  charge/discharge  mode.  The  size 
of  a  finite  volume  is  about  1/20  of  an  electrode  thickness.  A 
typical  discharge  simulation  requires  about  5  min  of  CPU 
time  on  a  PC  with  450  MHz  CPU. 
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4.  Experimental 

All  experiments  are  performed  on  a  commercial  nickel- 
metal  hydride  battery  of  nominal  capacity  95  Ah,  which 
consists  of  ten  cells  in  series.  Experiments  were  performed 
with  a  computer-controlled  battery  testing  system  (Arbin 
Instruments).  To  keep  the  isothermal  conditions,  fans  have 
been  used  to  induce  forced  convection,  and  the  variations  of 
temperature  were  determined  to  be  within  5  °C  by  thermo¬ 
couples  placed  on  the  surface  of  the  tested  batteries.  The 
same  charge  process  was  used  for  all  discharge  experiments: 
20  A  CC  step  until  the  battery  voltage  reaches  14.35  V, 
followed  by  a  CV  step  until  the  charge  capacity  reaches 
100  Ah,  finally  rests  for  1  h  after  each  full  charge.  The 
discharge  cutoff  voltage  for  all  low  rates  is  10.5  V,  but  this 
is  lowered  to  10.0,  9.5,  and  8.0  V  for  1,  2,  and  3  C-rates, 
respectively.  For  all  charge  experiments,  fully  discharge  the 
battery  by  constant  0.1  C-rate  until  the  cutoff  voltage  of 
10.5  V.  During  charge,  the  cutoff  voltage  was  set  at  15  Vat 
high  rates,  and  at  low  rates  batteries  were  allowed  to  reach 
the  maximum  voltage. 


5.  Results  and  discussion 

5.1.  Experiment  data 

5.1.1.  Reproducibility  of  battery  data 

The  mathematical  model  is  developed  to  predict  and 
simulate  the  behavior  of  a  single  cell.  Therefore,  only  when 
the  variations  among  the  cells  in  a  battery  module  were  small 
enough,  the  simulation  could  be  used  to  predict  the  battery 
performance  and  direct  its  design.  Fig.  2  shows  the  constant 
current  discharge  curve  and  error  bars  for  different  cells  at 
0.2  C-rate  discharge.  It  is  obvious  that  the  performances  of 
different  cells  are  quite  similar.  Therefore,  it  is  reasonable  to 


use  the  cell  model,  which  assume  that  each  cell  has  the  same 
behavior  as  others,  to  simulate  an  entire  battery  and  predict 
its  behaviors. 


5.1.2.  Determination  of  the  surface  area  of  electrode 
The  surface  areas  of  both  electrodes  are  determined  by 
BET  method;  for  the  Ni  electrode,  it  has  a  BET  surface  area 
of  26.5  m2/g,  and  for  the  MH  electrode,  it  has  a  surface  area 
of  1.9  m2/g.  The  specific  surface  area  of  each  electrode,  a0, 
can  be  expressed  by: 

a0  =  pS  (9) 

where  S  is  the  BET  surface  area  and  p  the  density  of  each 
electrode.  Assuming  that  all  particles  are  spherical,  thus  the 
total  solid  volume  is  given  by: 

Vs  =  V(1  —  e0)  =  |«7ir3  (10) 

The  total  solid  surface  area,  s  is  expressed  by: 
s  =  Amir2  (11) 


Therefore,  the  surface  area  a0  can  be  related  to  porosity  and 
the  radius  of  particles  by  the  following  equation: 


3(1 -eo) 

r 


(12) 


Assume  a  known  porosity,  the  equivalent  radius  can  be 
determined  from  the  above  equation.  All  geometrical  para¬ 
meters  that  are  different  from  those  in  Gu  et  al.  [12]  are  now 
listed  in  Table  1. 


5.1.3.  Oxygen  evolution 

It  is  known  that  the  effect  of  oxygen  reaction  is 
important  at  the  late  stage  of  charge.  In  order  to  describe 
the  oxygen  reaction  accurately,  constant  potential  experi¬ 
ment  is  carried  out,  in  which  the  battery  is  first  charged  at 
C/4  rate  until  a  sufficiently  high  potential  to  ensure  nearly 
full  charge,  and  subsequently  the  potential  is  held  steady 


Table  1 


Cell  geometrical  parameters 


Parameter 

Value 

Unit 

Ni  electrode 

Thickness 

0.045 

cm 

Projected  electrode  area 

14  x  11 

cm2 

rNi 

784000 

cm2/cm3 

^Ni 

0.33 

flNi 

1.56  x  10~4 

cm 

Separator 

Thickness 

0.017 

cm 

MH  electrode 

Thickness 

0.029 

cm 

aMH 

57490 

cm2/cm3 

£mh 

0.52 

rMH 

2.5  x  10~5 

cm 
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until  the  charge  current  reaches  a  steady  value.  This 
current  can  be  considered  as  the  oxygen  reaction  current. 
Due  to  the  small  transfer  current  density,  it  is  suitable  to 
assume  that  all  concentration  distributions  within  the 
battery  are  uniform,  and  that  the  ohmic  drop  as  compared 
to  the  kinetic  polarization  within  the  porous  electrode  is 
small,  namely,  the  potential  distribution  within  the  elec¬ 
trode  is  uniform.  Based  on  these  assumptions,  Tafel 
equation  can  be  used  to  describe  oxygen  evolution  reac¬ 
tion  as  shown  below: 


t]  =  a  +  b  log(i)  (13) 

where  1 7  is  the  overpotential  of  the  oxygen  evolution  reaction 
and  i  the  true  reaction  current  density  of  the  oxygen  reaction, 
which  can  be  calculated  by: 

*  =  \  (14) 

ciqL 


Here  is  is  the  superficial  current  density,  i.e. 
I 


(15) 


l  is  the  thickness  of  the  Ni  electrode,  I  is  the  total  current 
passing  through  the  battery,  and  A  is  the  projected  surface 
area  of  the  Ni  electrode.  The  values  of  a  and  b  are  found  in 
Fig.  3.  If  the  ohmic  drop  cannot  be  neglected  or  the  con¬ 
centration  distribution  is  not  uniform,  porous  electrode 
theory  [14]  should  be  used  instead.  From  the  following 
two  equations: 


2.3  RT 

a-  „  log  '0 

(16) 

aaF 


2.3  RT 
a  aF 


(17) 


the  exchange  current  density  of  oxygen  evolution  reaction, 
io  =  7.83  x  1 0  1  s  A/cm2,  and  the  transfer  coefficient  of 
anode  reaction,  y.a  =  1.5,  are  found. 


Fig.  4.  Comparison  of  different  equilibrium  potentials. 


5.2.  Model-experimental  comparisons 


The  difference  between  the  predicted  cell  open-circuit 
potential  and  experimental  data  is  clearly  seen  in  Fig.  4, 
where  the  low  rate  (0.05  C-rate)  data  is  compared  to  various 
thermodynamic  expressions.  Note  that  Fig.  4  assumes  that  a 
low  rate  experiment  can  be  used  to  approximate  the  equili¬ 
brium  potential  of  the  cell  as  a  function  of  SOC.  In  other 
words,  it  is  assumed  that  kinetic,  ohmic,  mass  transfer  losses 
and  self-discharge  due  to  oxygen  evolution  can  be  neglected 
under  such  rate. 

In  the  Ni-MH  battery  system,  the  capacity  of  MH  elec¬ 
trode  is  much  larger  than  that  of  Ni  electrode  due  to  the 
overcharge  and  overdischarge  reserve.  Therefore,  the  state 
of  charge  of  the  MH  potential  usually  stays  at  the  middle 
part.  According  to  Haran  et  al.’s  [8]  experimental  data,  the 
equilibrium  potential  of  MH  electrode  in  this  part  is  nearly 
constant.  Thus,  during  discharge/charge  the  MH  potential 
can  be  assumed  constant  at  low  rates  (0.05  C-rate  for  dis¬ 
charge  and  0. 1  C-rate  for  charge).  Accounting  for  the  inter¬ 
calation  in  the  Ni  electrode,  a  one-parameter  Margules 
model  was  used  in  Gu  et  al.  [12]  to  correct  the  activity  of 
proton.  A  two-parameter  Margules  expression  correcting 
Nernst  equation  for  Ni  electrode  is  also  shown  in  Fig.  4, 
which  can  be  expressed  as: 


E° 


0.3854  -) - 1 — Inf— — ^ 

38.94  VI  -x) 

I9W  (2  x0'8 186— 0-8186  x  4(1  -  x)  -2  x  0.9498 
+  6  x  0.9498(1  -  x)  -  3  x  0.9498(1  -  x)2)  (18) 


where  l*'  is  the  equilibrium  potential  of  the  Ni  electrode  and 
x  is  the  state  of  charge.  It  can  be  seen  that  the  two-parameter 
Margules  model  makes  the  curve  flatter  than  the  one  pre¬ 
dicted  by  Nernst  equation.  However,  the  slope  of  the  curve  at 
intermediate  SOC  is  much  steeper  than  that  seen  experi¬ 
mentally,  and  it  still  cannot  capture  the  main  shape  of  the 
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equilibrium  potential  curve.  Considering  the  inability  of 
these  theoretical  models  to  fit  the  data,  the  data  was 
curve-fitted  to  the  following  equation: 

E°  =  0.314  +  — —  x  In  ( — — — ) 

38.94  VI  -x) 

+  0.5997x-  1.1889x2  +  0.5937x3  +  0.0985-y4  (19) 

This  expression  is  subsequently  employed  in  the  mathema¬ 
tical  model. 

According  to  Paxton  and  Newman  [10]  and  Gu  et  al.’s 
[12]  work,  the  proton  diffusion  would  not  be  important  if  the 
proton  diffusion  length  is  small.  The  length  measured  in  this 
work  is  substantially  smaller  than  the  length  used  by  Gu  et  al. 
[12].  Therefore,  the  present  diffusion  coefficient  of  proton  is 
the  same  as  Gu  et  al.  [12],  measured  by  MacArthur  [15]. 

From  Paxton  and  Newman’s  [10]  discussion,  the  para¬ 
meters  of  the  nickel  electrode  affect  the  battery  discharge 
behavior  most.  Thus,  we  fit  the  highest  rate  discharge  curve 
to  obtain  a  more  realistic  set  of  kinetic  parameters  of  Ni 
electrode.  Table  2  shows  new  kinetic  parameters  used  in  the 
corresponding  simulation  study.  Note  that  the  product  of  the 
specific  surface  area  and  exchange  current  density  used  in 
the  present  work  is  only  slightly  different  from  the  value 
used  by  Gu  et  al.  [12]  for  the  Ni  electrode  and  is  exactly  the 
same  for  the  MH  electrode.  In  addition,  the  transfer  coeffi¬ 
cient,  a,  which  affects  the  shape  of  the  discharge  curve  of  Ni 
electrode,  is  another  parameter  adjusted  for  better  agreement 
with  the  experimental  data. 

Fig.  5  displays  a  comparison  between  experimental  and 
simulated  discharge  curves  over  a  wide  range  of  C-rates.  It 
can  be  seen  that  the  main  shapes  of  the  experimental  dis¬ 
charge  curves  were  captured  successfully  by  the  simulation. 
Especially  at  0.1  C-rate,  even  the  initial  quick  drop  in  the  cell 
potential  was  predicted  well,  showing  that  at  low  rate  dis¬ 
charge  Nemst  equation  controls  the  behavior  of  the  battery.  At 
the  middle  part,  all  predictions  are  in  good  agreement  with  the 
measured  curves.  It  can  be  found  that  the  prediction  of  battery 
capacity  is  fairly  accurate  at  all  rates.  Changing  the  diffusion 
coefficient  in  the  model  does  not  affect  the  discharge  curve, 
thus  indicating  that  the  proton  transport  is  not  a  limiting  step, 
in  accordance  with  the  observations  of  Paxton  and  Newman 
[10]  and  Gu  et  al.  [12].  In  summary,  the  present  Ni-MH 
model  coupled  with  the  well-characterized  kinetic  and  trans¬ 
port  properties  results  in  generally  satisfactory  agreement 
with  discharge  data  over  a  wide  range  of  rates. 

Due  to  the  hysteresis  effect  on  Ni  electrode,  the  open- 
circuit  potential  equation  for  both  discharge  and  charge 


Table  2 

Kinetic  parameters  for  Ni  electrode 


Parameter 

Present  paper 

Gu  et  al.  [12] 

Unit 

$ 

0.73 

0.5 

a? 

0.27 

0.5 

“‘o‘ 

0.4666 

0.4194a 

A/cm3 

a  Obtained  from  individual  values  of  a  and  /[[  listed  in  Gu  et  al.  [12], 


Fig.  5.  Comparison  of  simulated  battery  voltages  (solid  lines)  with 
experimental  data  (symbols). 

processes  may  not  be  the  same.  Because  of  oxygen  evolution 
reaction  during  the  charge  process,  which  can  significantly 
affect  the  end  part  of  the  charge  process,  we  cannot  use  the 
same  method  as  for  discharge  to  determine  the  open-circuit 
equilibrium  potential  of  the  charge  process.  Instead,  we  use 
the  same  formula  as  the  discharge  part,  but  slightly  adjust  the 
coefficients  to  obtain  the  charge  open-circuit  potential  of  Ni 
electrode  as  follows  (also  shown  in  Fig.  6): 

£°  =  0.494  +  — —  x  Inf -3— ) 

38.94  Vl  -x) 

+  0.0784b:  -  0 . 224424.x2 + 0 . 044965.*3  +  0.044965x4 

(20) 

With  the  above  open-circuit  potential  for  charge  but  otherwise 
all  the  same  physico-chemical  properties  as  for  discharge, 


Fig.  6.  Open-circuit  potential  of  Ni  electrode  with  respect  to  state  of 
charge  (solid  lines  for  model  and  symbols  for  experimental  data). 
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Fig.  7.  Comparison  of  simulation  battery  voltages  (solid  lines)  with 
experimental  data  (symbols). 

Fig.  7  displays  a  comparison  of  the  simulated  and  measured 
charge  curves  over  a  range  of  C-rates.  It  can  be  seen  that  the 
simulations  agree  well  with  the  experimental  data.  From  the 
plot,  we  can  see  that  at  the  beginning  of  charge,  the 
potentials  for  all  rates  increase  fast,  then  become  nearly 
constant  until  90  Ah,  followed  by  a  steep  increase  again  near 
full  charge,  and  lastly  a  constant  plateau  due  to  the  balance 
of  oxygen  evolution  and  recombination  reactions  within  the 
battery.  It  is  interesting  to  observe  that  the  plateau  of  each 
curve  in  overcharge  occurs  almost  at  the  same  charge 
capacity.  This  seems  to  suggest  that  the  oxygen  cycle  does 
not  start  until  full  charge. 

From  the  model,  one  can  conveniently  calculate  the 
respective  distributions  of  the  main  reaction  and  oxygen 
evolution  reaction  currents.  It  can  be  seen  from  Fig.  8  that 
the  side  reaction  current  is  relatively  small  before  the  charge 


Capacity  (Ah) 

Fig.  9.  Comparison  of  battery  charge  curves  with  (solid  lines)  and  without 
oxygen  reaction  (dashed  lines). 

capacity  reaches  85  Ah  for  all  rates  charge.  After  that  the 
side  reaction  rate  increases  sharply  until  reaching  the  total 
charge  current.  Charging  beyond  102  Ah,  nearly  the  entire 
charge  current  for  all  rates  is  to  supply  the  oxygen  reaction. 

Fig.  9  shows  the  predicted  effects  of  oxygen  reaction  on 
the  charge  potential.  During  charge,  most  part  of  the  charge 
curve  is  not  affected  by  oxygen  reaction,  in  consistency  with 
Fig.  8.  It  is  interesting  that  oxygen  side  reaction  affects 
charge  potential  at  low  rates  much  more  than  at  high  rates. 
This  is  the  main  reason  why  it  is  not  suitable  to  consider  the 
low  rate  charge  curve  as  the  equilibrium  potential  due  to  the 
oxygen  evolution. 

5.3.  Battery  design 


Fig.  8.  Comparison  of  oxygen  reaction  current  (dashed  lines)  and  main 
reaction  current  (solid  lines)  for  various  C-rates  during  charge. 


The  charge  algorithm  is  an  important  part  of  battery 
applications,  especially  for  electric  and  hybrid  electric 
vehicles,  because  it  is  desirable  to  charge  the  battery  in 
the  shortest  time  and  at  the  highest  charge  efficiency.  How¬ 
ever,  both  goals  cannot  be  achieved  simultaneously  because 
the  parasitic  current  associated  with  oxygen  reaction 
increases  with  charge  rate  and  state  of  charge.  Charge 
algorithms  are  typically  optimized  by  experimental  trial- 
and-error.  Using  a  validated  model  to  assist  in  this  optimiza¬ 
tion  can  be  a  cost-effective  alternative. 

Charge  algorithms  typically  consist  either  of  a  constant 
current  step,  or  of  a  constant  current  step  followed  by  a 
constant  voltage  step.  Fig.  10  displays  the  latter  method  for 
the  12  V  Ni-MH  battery  and  shows  the  final  state  of  charge 
and  charge  efficiency  as  functions  of  the  voltage  step  level. 
The  charge  efficiencies  for  all  charge  rates  decrease  with  the 
final  voltage  level,  due  to  the  side  reactions  having  more 
effect  at  higher  potentials.  The  final  SOC  is  solely  deter¬ 
mined  by  the  voltage  level,  regardless  of  the  rate  of  current 
used  in  the  first  step,  and  increases  with  the  voltage  level. 
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Fig.  10.  Charge  efficiency  (solid  lines)  and  final  state  of  charge  (dashed  lines)  for  a  12  V  Ni-MH  battery  by  constant  current  followed  by  constant  voltage 
method. 


The  charge  efficiency  shows  the  opposite  trend;  the  effi¬ 
ciency  decreases  at  a  certain  C-rate  as  the  voltage  level  is 
increased.  On  the  other  hand,  the  figure  shows  that  the 
charge  efficiency  increases  with  charge  rate  at  a  certain 
charge  voltage.  This  increase  is  shown  more  clearly  in 
Fig.  11,  which  further  indicates  that  the  charge  efficiency 
becomes  nearly  constant  above  0.5  C.  It  is  evident  from  the 
discussion  above  that  low  rate  charging,  e.g.  0.1  C,  up  to  a 
particular  constant  voltage  level  is  not  desirable  due  to  the 
low  charge  efficiency.  A  compromise  between  these  com¬ 
peting  effects  must  be  made  to  determine  the  optimum 
algorithm. 


Fig.  11.  Charge  efficiency  as  a  function  of  C-rate  by  constant  current 
followed  by  constant  voltage  method. 


Fig.  12  shows  the  charge  efficiency  with  respect  to  the 
charge  input  for  constant  current  charge  method.  The  charge 
efficiency  increases  with  the  charge  rate  and  then  reaches  a 
maximum  between  0.5  and  1  C.  The  efficiency  decreases 
with  charge  rates  above  1  C.  The  charge  efficiency  is  high  at 
the  beginning  of  charge  for  low  charge  rates  (e.g.  <0.25  C), 
but  then  decreases  quickly  with  the  charge  input.  It  is 
interesting  that  the  charge  efficiency  at  low  charge  rates 
for  constant  current  charging  is  lower  than  for  higher  rates. 
These  observations  show  that  charge  rates  near  1  C  are 
optimal  for  achieving  maximum  charge  efficiency. 
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Fig.  13  displays  the  oxygen  reaction  current  with  respect 
to  charge  input  for  different  charge  rates.  As  expected,  larger 
charge  rates  result  in  higher  side  reaction  currents.  This 
result  seems  to  suggest  that  the  charge  efficiency  decreases 
with  the  increase  of  charge  rate  due  to  the  larger  side 
reaction  current  at  the  higher  charge  rate.  However,  the 
charge  time  is  also  different  for  different  rates,  and  the  total 
charge  capacity  that  goes  into  the  side  reaction  is  determined 
by  the  product  of  time  and  current.  The  charge  capacity 
consumed  by  the  side  reaction  at  different  charge  rates  with 
respect  to  charge  input  is  plotted  in  Fig.  14.  The  side  reaction 
charge  capacity  is  the  lowest  at  intermediate  charge  rates, 
with  the  highest  for  both  0. 1  and  4  C,  which  agrees  with  the 
conclusion  drawn  from  Fig.  12. 


Fig.  14.  Charge  capacity  consumed  by  the  oxygen  reaction  for  different 
charge  rates. 


Fig.  15.  The  final  state  of  charge  as  a  function  of  the  terminating  voltage 
in  constant  current  charge. 

Fig.  15  shows  the  relationship  between  the  final  state  of 
charge  and  the  terminating  battery  voltage  during  constant 
current  charging.  The  relationship  for  all  charge  rates  exhi¬ 
bits  the  same  trend:  initially  the  state  of  charge  increases 
gradually,  followed  by  a  sharp  increase,  and  then  tails  off  as 
the  SOC  approaches  100%.  The  charge  efficiency  for  a  given 
charge  rate,  however,  decreases  with  the  terminating  voltage 
as  shown  in  Fig.  16.  The  main  reaction  charge  efficiency 
drops  off  sharply  as  the  cell  potential  approaches  a  certain 
value  because  the  charge  current  is  used  dominantly  in  the 
oxygen  cycle.  Figs.  15  and  16  can  be  used  to  determine  the 
cutoff  voltage  of  constant  current  charge  algorithm  accord¬ 
ing  to  the  desired  charge  efficiency  and  state  of  charge. 
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6.  Conclusion 

An  approach  to  closely  coupling  modeling  with  experi¬ 
mentation  has  been  attempted  for  a  commercially  available 
Ni-MH  battery.  The  battery  parameters,  the  equilibrium 
potential  of  the  Ni  electrode  and  its  hysteresis  effect,  and 
the  kinetics  of  oxygen  reaction  were  characterized  and  incor¬ 
porated  into  a  previously  developed  first-principles  model  to 
result  in  a  generally  good  agreement  between  experimental 
and  predicted  discharge  as  well  as  charge  curves.  The  vali¬ 
dated  mathematical  model  was  then  used  to  study  charge 
algorithms.  Low  rate  charging  is  not  desirable  due  to  the 
lower  charge  efficiency  for  both  charge  algorithms:  the  con¬ 
stant  current  step  and  the  constant  current  step  followed  by 
constant  voltage  step.  The  model  predictions  show  that  1  C- 
rate  charge  can  be  used  to  quickly  charge  the  battery,  while 
still  achieving  the  highest  charge  efficiency.  Future  work  will 
concentrate  on  the  thermal  effects  in  Ni-MH  batteries. 
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